Singular continuous spectra in a pseudo-integrable billiard 
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The pseudo-integrable barrier billiard invented by Hannay and McCraw [J. Phys. A 23, 887 
(1990)] - rectangular billiard with line-segment barrier placed on a symmetry axis - is generalized. 
It is proven that the flow on invariant surfaces of genus two exhibits a singular continuous spectral 
component. 
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PACS number: 05.45.-a 

Billiards have attracted much attention in the recent 
decades as simple dynamical systems in classical and 
quantum mechanics. The free motion of a particle in 
a domain within a hard elastic boundary shows a large 
spectrum of behaviour depending on the shape of the 
boundary, ranging from complete regularity to various 
degrees of chaoticity. 

Polygonal billiards are neither chaotic the 
Kolmogorov-Sinai entropy is zero - nor integrable in 
the sense of Liouville-Arnol'd Jl| (apart from the rectan- 
gles, the equilateral triangles, the 7r/2, 7t/4, 7r/4-triangles 
and the tt/2, it/3, 7r/6-triangles). The motion inside a 
typical polygon is conjectured to be ergodic on the 
three-dimensional constant-energy surfaces in phase 
space, while the motion inside a rational polygon (all 
angles are rationally related to it) is restricted to two- 
dimensional invariant surfaces, like in integrable systems, 
but the genus of the surfaces is larger than 1. Rational 
polygonal billiards are therefore characterized as pseudo- 
integrable P). It is rigorously proven that the flow on 
these surfaces is ergodic and not mixing |3|. Ergodicity 
means that a particle with a typical initial momentum 
explores the entire invariant surface (which is the bil- 
liard table in the configuration space projection). Since 
the mixing property is excluded, a small cluster of par- 
ticles with identical initial momenta cannot spread out 
uniformly in time on the invariant surface. But a nonuni- 
form spreading which allows occasional reclustering with 
decreasing frequency - weak mixing in mathematical 
terms (see e.g. g) — is possible. In fact, weak mix- 
ing is believed to be generic, but this is only shown for 
polygons all of whose sides are horizontal or vertical. Nu- 
merical evidence for weak mixing has been reported for 
the "square-ring billiard" ||; 

see also related numerical 
studies on rational || and irrational triangles . 

Weak mixing (in the absence of the stronger mix- 
ing property) implies interesting spectral properties 
More precisely, the Fourier transform of a typical func- 
tion on the phase space with respect to time has a fractal- 
like structure in the frequency domain. Such a spec- 
trum is characterized as singular continuous, in contrast 
to discrete spectra of integrable systems (nonmixing, 
quasiperiodic motion) and absolutely continuous spectra 
of chaotic systems (mixing motion). Singular continu- 



ous spectra are well known in other fields of physics to 
appear at the border between integrability and chaos, 
e.g. strange nonchaotic attractors (see, e.g., fl0| , |il"ll ) , the 
(kicked) Harper model Jh2[ and models for nonperiodic 
atomic structures (see |13| and references therein). 

The spectral properties of a phase space function are 
related to its autocorrelation function (AF). For discrete 
spectra this function is (quasi-) periodic. In the periodic 
case the (normalized) AF returns exactly to 1, while in 
the quasiperiodic case it comes arbitrarily close to 1. An 
exponentially decreasing AF (for increasing time) indi- 
cates an absolutely continuous spectrum. An AF cor- 
responding to a singular continuous spectrum does not 
usually decay to zero, but also does not return to 1. 

A nongeneric subclass of pseudo-integrable systems 
is given by almost-integrable billiards G3j. A member 
of this class is, roughly speaking, composed of several 
copies of a single completely integrable billiard, e.g. the 
7r/3, 27r/3-rhombus Eq] consists of two identical equilat- 
eral triangles. Motion on invariant surfaces of almost- 
integrable billiards is not weakly mixing . The aim of 
this letter is to establish that the spectra of these systems 
nevertheless can have a singular continuous component 
(together with a discrete component). We prove this for 
a simple model, a generalized version of the "barrier bil- 
liard" [[l6| . It is the first time that the existence of singu- 
lar spectra in a given billiard is rigorously shown. This is 
achieved by decomposing the dynamics into a continuous 
and a discrete part: the quasiperiodical motion in the in- 
tegrable subbilliard and the transitions between different 
copies. The latter type of motion, responsible for the sin- 
gular continuous component, is discussed with the help 
of the AF. 

We consider a point particle with unit mass elastically 
bouncing inside a polygonal enclosure consisting of a ver- 
tical line of length 1 — j3 € (0, 1) symmetrically placed in 
a rectangle with width L G (0, oo) and normalized height 
1, see Fig. |l](a). The special case (3=1/2 corresponds to 
the usual barrier billiard [M. The trivial energy depen- 
dence of the system is ignored by setting the energy to 
1/2, or equally the magnitude of the momentum (p x ,Py) 
to 1. Given an initial momentum only a finite number 
of directions is realized during the time evolution, due to 
the fact that all angles in the polygon are rational mul- 
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tiples of 7r. In other words, the motion in phase space 
takes place on two-dimensional invariant surfaces. Us- 
ing the general formula for the genus of such surfaces ||] 
gives 2, i.e. the surfaces have the topology of two-handled 
spheres and not that of tori (single-handled spheres). 
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FIG. 1. (a) Barrier billiard, rectangle with a vertical line con- 
necting the origin of the coordinate system (x, y) = (0, 0) with the 
point (x, y) = (0,1-/3). (b) Symmetry reduced version. 

Our model satisfies the definition of almost- 
integrability given by Gutkin (m) for all values of L 
and (3. It is composed of two copies of an integrable 
rectangular billiard, which can be obtained by sym- 
metry reduction, see Fig. 0(b). The dynamics can be 
split into two parts. First, the quasiperiodical dynam- 
ics (\x(t)\,y(t)) on the invariant tori of the symmetry 
reduced system. Second, the dynamics of the sign of x, 
s{t) = ±1. 

The quasiperiodical component is most elegantly de- 
scribed in action-angle variables. The actions I x — 
\p x \L /2tt and I y — \p y \/ir label the invariant tori while 
the angles (divided by 2n and computed modulo 1) de- 
scribe the flow on a given torus 
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tt%. The 



with the frequencies u> x = 4ir 2 I x /L 2 and u> y 
flow on the torus is ergodic if and only if the winding 
number p = uj v /uj x = I y L 2 / (AI X ) is irrational. The mo- 
tion in configuration space is explicitly given by 



y(t) = f(4> y (t) 



1-/3 



(mod 1)) 



(3) 
(4) 



with the piece-wise linear function 



/(«) = 



2u 



if 0<u<l/2 



2(1 — u) otherwise. 



The shift in the argument of the function / in Eq. (^) 
will be convenient later. 



Now we consider the non-quasi-periodical part of the 
dynamics. Clearly, the value of s(t) changes when the 
trajectory impinges on the segment {x — and 1— (3 < y} 
(or {(f) x = and < <f) y < /3}). We therefore introduce 
the Poincare section 4> x (nT) — with T — 2tt/uj x and 
the integer n. We obtain 



b y ,„+i = (j> v , n + p (mod 1) 

S n +l = Sn$(4>y,n) 



(5) 
(0) 



with 



$(0) 



-1 if < <f) < (3 
1 otherwise. 



Note that diffractive orbits, i.e. orbits which hit the edge 
of the barrier, can be neglected since they have measure 
zero in phase space. The two-dimensional mapping (|J)- 
(p) has the solution 



s n = s o n™ =o $(0 a! o + jp (mod 1)). 



(7) 



This equation together with Eq. (|l])-(Q) and n given by 
the integer part oit/T solve the billiard problem. But the 
solution is only formal, because in order to compute (|^) 
one has to iterate (||)-(||). A similar, but more com- 
plicated, solution of the rhombus billiard was reported 
in Jl5[ . But its formality was not recognized since the au- 
thors focussed only on the equation corresponding to (||) , 
while the equation corresponding to (|^) was hidden in the 
notation. 

Having defined the dynamical system, we now dis- 
cuss the spectral properties of a typical smooth 
phase space function F(x,p x , y,p y ) (nontypical functions 
F(x 2 ,p 2 ,y,p y ) obviously have discrete spectra). We in- 
vestigate firstly the function F(t) = x(t) (the argumen- 
tation for p x (t) is identical). We will see that x(t) has a 
singular continuous spectrum, which implies that a typ- 
ical function has a mixture of discrete and singular con- 
tinuous spectrum. With x(t) = s(t)|a;(t)| it is clear that 
the Fourier transform of x(t), 



x(t)e 



(8) 



can be written as convolution of the Fourier transforms 
of \x{t)\ and s(t), X{lu) and S{u) 



(9) 



Since \x(t)\ is periodic with period T we have 

oo 

X(ui') = ^2 A m 5{uj' -muj x ). 



m— — oo 



Substituting into Eq. (M) gives 
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muj x ) 



(10) 



m— — oo 



Since s[t) changes in discrete steps it is suitable to 
consider the discrete-time Fourier transform ol the se- 
quence [s n ] 



which is related to the continuous-time Fourier transform 
via 



S(u) = — (l-e 4wT )E(u;T). 
iuj 



(11) 



Note that the prefactor goes to the constant T as ui — > 0. 
We can now take advantage of the fact that the map- 
ping (H)-(U) is well studied. It is proven (see |nj and ref- 
erences therein) that if p is irrational then S(w) is singu- 
lar continuous for almost all (3 (the spectrum is discrete, 
e.g., if = np (mod 1)). Since this holds trivially also 
for any nonconstant function of s(t) we conclude that the 
s-motion is weakly mixing The fact that almost all 
trajectories in our billiard have irrational winding num- 
ber p, independent of the values of the parameters L 
and /?, ensures that typical trajectories give rise to sin- 
gular spectra S(w). Together with Eq. ( |l0| ) and ( |Tl| ) this 
implies that X{lo) is typically singular continuous. Nev- 
ertheless, there exist nonconstant functions of x(t), e.g. 
x 2 (i), which do not have a singular continuous spectrum, 
so the x-motion cannot be weakly mixing. This result in- 
dicates that a purely numerical detection of weak mixing 
by investigating the spectral properties of just a single 
phase space variable may lead to wrong conclusions. 

Figure || shows an example of |X(w)| 2 close to the fre- 
quency lu x (so essentially \S(ui — u) x )\ is shown) for the 
case (3 — 1/2 and p equal the reciprocal of the golden 
mean, p gm = (\/5 — l)/2. The integration (^) has been 
performed with a fading factor exp ( — 2i/t max ) and trun- 
cated at t max = 10000. The enlargement in Fig. || re- 
veals some selfsimilarity of S(lu) for this set of parame- 
ters. The minor lack of symmetry is due to the prefactor 
in Eq. (|TT|); S(lj) is perfectly selfsimilar as shown by 
Feudel et al. [Tl| ] by using a renormalization group ap- 
proach (compare the magnification in Fig. |^ with Fig. 3 
in ]ll[). Feudel et al. have also calculated the general- 
ized dimensions: (i) the Hausdorff dimension D = 1, i.e. 
the spectrum is continuous, (ii) The information dimen- 
sion D\ « 0.75, i.e. the spectrum is singular, (hi) The 
correlation dimension Di ~ 0.65, which is the exponent 
of the power-law decay t~ D2 of the integrated AF Q. 
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FIG. 2. Power spectrum |Jf(a))| 2 for f3 = 1/2 and p = Pgm- 
The dashed line marks the fundamental frequency of |x(i)|, 
u) x ?a 3.95187. The inset shows a magnification. 

Before discussing the integrated AF, let us first con- 
sider the (normalized) AF (x(t) has zero mean value) 



R x (t) = 



(x(t + r)x(t)) 



(12) 



where the average (. . .) is taken over t. We restrict our- 
self to the times r = nT . By using the periodicity of 
|a:(t)| it is an easy matter to show that R x (nT) is equal 
to R s (n), the AF of the process s n (where the average is 
over the discrete time n). This means we discuss here the 
correlations of the transitions between the two copies of 
the integrable subbilliard. For the special case (3 = 1/2, 
R s (n) as function of n is studied in Jl(| for different p. It 
is shown that R s (n) vanishes for odd n. For large even 
n it is mostly close to zero but has large peaks at posi- 
tions which roughly grow exponentially. For example, if 
P — Pgm then peaks with amplitude « 0.55 are located at 
all even Fibonacci numbers surrounded by smaller peaks 
in a selfsimilar manner. 

We have determined the integrated AF 
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(13) 



up to n max = 10000. Its asymptotic behaviour is found 
to be in agreement with the expected power-law decay, 
which allows us to compute Di numerically. For fixed 
P = Pgm, D2 as function of (3 (using 1000 randomly 
distributed /3-values) has mean value 0.458 and shows 
very strong fluctuations, ranging from large decay rates 
Z?2 = 0.68 to zero decay. These fluctuations are by no 
means surprising, since there are dense sets of /3-values 
(e.g. (3 = np (mod 1)) for which D2 vanishes. We have 
also calculated D2 for fixed (3 = 1/2 as a function of 
p (using 1000 p- values). Here, D2 is more or less uni- 
formly distributed between zero and 0.94 around a mean 
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value of 0.469. Due to the sensitive dependence of D2 on 
the winding number, it seems useful to perform an aver- 
age of the AF over all invariant surfaces (i.e. over all p) 
in order to characterize the correlation properties of the 
transitions between the two copies of the subbilliard as a 
whole. But this kind of averaging does not yield a reason- 
able correlation function since for an integrable system 
such a function would indicate correlation decay p9[ . 

The fact that the integrated AF (|l^) decays to zero 
but R s does not confirms the weak-mixing property of 
the s-motion. Because of R x (nT) = R s (n) one could 
naively draw the wrong conclusion that the x-motion is 
also weakly mixing. This confusion can be resolved if 
one clearly distinguishes between continuous-time and 
discrete-time approaches. Trivially, R x does not de- 
cay for large continuous times r if it does not decay 
for large discrete times n. But we expect that the in- 
tegrated (now really an integration not a summation) 
AF of x(t) does not decay, in contrast to the integrated 
AF of the process s„ (the integrated AF of a general 
function F(x, y,p x ,Py) obviously shows zero decay rates, 
since y(t) is periodic). This agrees with results of the 
study on an almost-integrable version of the "square-ring 
billiard" Q. It is worth emphasising that the discrete- 
time approach, which is always used in numerical studies 
where integrations are replaced by summations (for bil- 
liards the discrete-time approach is often used explicitly) , 
can wrongly indicate weak mixing, even though this is 
quite unlikely. 

In terms of the billiard motion we have the following 
scenario. The motion is not mixing and even not weakly 
mixing; that means, a small cluster of particles with the 
same initial momenta will not spread out on the entire 
billiard table. But the motion is "partly weak mixing" 
in the sense that after some transient time (when the AF 
is close to zero) there are two clusters, one in the left 
and one in the right part, each of them consisting of ap- 
proximately half of the particles. At "resonant" times 
(corresponding to the peaks of the AF) one of the two 
clusters contains more particles than the other one. But 
these reclusterings become more and more rare with in- 
creasing time. 

To summarize, we have proven that the motion in the 
pseudo-integrable barrier billiard exhibits a discrete and 
a singular continuous spectral component. We expect 
that this is true for all almost-integrable billiards. 
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